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Abstract. We report on a detailed analysis of the cor- 
relation between the optical-UV and X-ray luminosities 
of quasars by means of Monte Carlo simulations, using a 
realistic luminosity function. We find, for a quasar pop- 
ulation with an intrinsically constant, mean X-ray loud- 
ness aox, that the simulated aox~~Lo relation can exhibit 
various 'apparent' properties, including an increasing aox 
with Lo, similar to what has been found from observa- 
tions. The determining factor for this behavior turns out 
to be the relative strength of the dispersions of the lu- 
minosities, i.e. their deviations from the mean spectral 
energy distribution at the optical and X-ray bands, such 
that a dispersion larger for the optical luminosity than for 
the X-ray luminosity tends to result in an apparent corre- 
lation. We suggest that the observed aox-Lo correlation 
can be attributed, at least to some extent, to such an ef- 
fect, and is thus not an underlying physical property. The 
consequences of taking into account the luminosity disper- 
sions in an analysis of the observed luminosity correlations 
is briefly discussed. We note that similar considerations 
might also apply for the Baldwin effect. 
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1. Introduction 

A study of the dependence of the spectral energy distri- 
bution (SED) of quasars on their luminosity and/or on 
cosmic epoch is particularly important for understanding 
the quasar phenomenon. In the optical-to-X-ray regime, 
the SED can be characterized by the broad band spec- 
tral index between 2500A and 2keV, which is defined as 
aox — —0.384 log(L2kcv/i °)- Quasars are known to 

2500A. 

exhibit strong luminosity evolution in the X-ray and the 
optical wave bands (e.g. Boyle 1994). However, there have 
been controversial discussions in the past as to whether 
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the evolution law is the same in these two energy bands. A 
dependence of aox on redshift or optical luminosity would 
indicate different evolution in the optical and the X-ray 
regime. Further, if aox depends on optical luminosity, this 
is equivalent to a non-linear relationship between X-ray 
and optical luminosity (Lx L"^, e 1). 

While most of the analyses agree on the result that 
aox is redshift independent, it has been claimed that aox 
increases with Lq, which implies that the objects with 
high optical luminosities are under-luminous in X-rays 
compared to their low luminosity counterparts (Avni & 
Tananbaum 1982, 1986; Kriss & Canizares 1985; Wilkes 
et al. 1994; Avni et al. 1995; Green et al. 1995). Generally, 
for a functional dependence of the form aox ~ PlogLo, a 
canonical slope of (3 0.1 was obtained, which is equiva- 
lent to a non-linear relation of the form Lx L^ "^ . 

Based on Monte Carlo simulations, Chanan (1983) 
suggested that a non-linear relation might arise even for 
an intrinsically linear dependence from observational flux 
limits and the large intrinsic scatter in the data. Thus, 
the observed aox - Lo correlation should not be considered 
as an underlying physical reality. The author also claimed 
that the choice of Lo as the independent variable is not jus- 
tified. However, Kriss & Canizares (1985) criticized these 
results by pointing out that they depend critically on the 
assumption of a Gaussian distribution for the luminosity 
functions. 

A study by La Franca et al. (1995) reinforced the idea 
of a linear relationship between the X-ray and the optical 
luminosity for quasars. They applied a regression algo- 
rithm to a large sample of quasars detected with Einstein, 
which accounts for errors in both variables and the intrin- 
sic scatter in the data, and found Lx on Lo- 
in a recent study of ROSAT detected quasars by 
Brinkmann et al. (1997), it has been shown by means of 
a simple argument that an apparent correlation between 
aox and log Lo can indeed emerge even for intrinsically un- 
correlated variables. Motivated by this idea, as well as by 
recent improvements concerning the shape of the quasar 
luminosity functions in the optical and the X-ray regime. 
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we carry out a detailed study of this controversial problem 
by means of a Monte Carlo analysis. We mostly use the 
logarithms of luminosities and denote them as 4c = log Lx 
and lo = logio- We use qo = 0.5, Hq = 50 kms~^Mpc~^ 
throughout this paper. All errors quoted are at the la 
level unless mentioned otherwise. 

2. Analysis of luminosity correlations 

2.1. Distribution of quasar luminosities 

We first describe the assumed parametric form of the mul- 
tivariate distribution '0(^0, 4c, z), which we use in the simu- 
lations. In previous studies (e.g. Avni & Tananbaum 1982, 
1986; Kriss & Canizares 1985) "0(^0, L, z) is commonly 
expressed as the product of the luminosity function (LF 
hereafter) of the primary luminosity (either 4c or lo) and 
the conditional distribution function of the secondary lu- 
minosity. The conditional distribution function determines 
the expected value of the secondary luminosity by assum- 
ing a functional dependence between the two luminosities 
and intrinsic dispersion. In previous models, the dispersion 
of aox was attributed to the dispersion in the secondary lu- 
minosity alone. In this paper we use a generalized descrip- 
tion of the multivariate distribution function. We assume 
that both luminosities show intrinsic dispersion instead of 
the secondary luminosity only. 

We introduce intrinsic optical and X-ray luminosities 
To and Ix- The observed luminosities lo and 4c are the in- 
trinsic luminosities modified by various mechanisms which 
produce a large scatter. The luminosities To and 7x are 
distributed according to their respective luminosity func- 
tions, (j)o(lo,z) and (j)y.{l^^z). Although lo and l^ are not 
directly observable, their ratio, i.e. the intrinsic aox, which 
is defined as 

aox = , ^° r = 0.384(Io - Ix), (1) 

can be approximated by the mean of the observed ckox 
distribution. As indicated by the observed correlation be- 
tween 4c and lo, Ix and lo are physically related. Assuming 
a redshift-independent relationship of the form 

Ix =/(Io)(xeIo , (2) 

we get an intrinsic dependence of cfox on lo, 

aox(Io) = Antio + const, where /3int = 0.384(1 - e). (3) 

If e = 1, we have /3int = 0, i.e. Sox is independent of the 
luminosities intrinsically. 

The observed optical and X-ray luminosities of a quasar 
are then 

^o=Io+<5/o, ^x=Ix + <5^x, (4) 

where Slo and Sl^ quantify the scatter around the intrin- 
sic luminosities in the optical and the X-ray band, respec- 
tively. The dispersion of c^ox can now be attributed to the 



dispersions in both 4c and lo. Assuming that Slo and Six 
are independent, we havea 

(Saoxf ^ 0.38A^{6lof + {Sl.n (5) 

where Saox is the dispersion of aox. 

Based on this generalized scenario, the multivariate 
distribution function ipi^o, Ix, z) can be replaced by a dis- 
tribution function depending on z. To, Ixj Slo and Sl^, which 
is given by 

1p*{lo,~lx,Slo,Sly,,z) = 

0o(Io, z) '5(Ix - fOo)) go{SloUo, z) gx{SlxUo, z), (6) 

where 4>o is the luminosity function for To, (5(7x — /(lo)) 
is a (5-function with / from Eq. ^ and go and are the 
conditional distribution functions for Slo and Six, respec- 
tively. In the following we describe these terms as they 
were implemented in the Monte Carlo analysis. 

We assume that both luminosity dispersions, Slo and 
51^, are independent of luminosity and redshift, and that 
go and g^ are given by Gaussian distributions with means 
of zero, and standard deviations Uo and a^, respectively. 
Thus, the distribution of aox is Gaussian with standard 
deviation 

aaox^0.38A{ao^ +ajy/^. (7) 

cFaox i® available from observational data, ranging from 
0.15 to 0.2 (e.g. Avni et al. 1995, Yuan et al. 1998). We 
define the parameter as the ratio of the standard devi- 
ations of the optical to the X-ray luminosity dispersions, 



Then, given and ctq/^^, cto and CTx can be determined 
using Eq.|^. 

Another input is the LF (j)^(lo,z), which is not di- 
rectly observable. Assuming a z-dependent power law for 
(j)o(lo,z) and a Gaussian distribution for go, it can be 
shown that the LF for the observed luminosity lo has the 
same functional form and evolution as that for To except 
close to the low-luminosity cutoff. Therefore, it is natu- 
ral to approximate (j)o(lo, z) by the observed optical lumi- 
nosity function (OLF) and its evolution. If Slo is small, 
(poQ-o, z) reduces to the observed OLF. 

We assume pure luminosity evolution for quasars. Us- 
ing the functional form of the OLF as given by Boyle 
(1994), which was derived from the UVX sample, we pa- 
rameterize 0o(Io, z) by a broken power law with 71 = — 1.6 
for Io< ^o and 72 = -3.9 for Io> 11, and Io(z) = 
Io(2=0)-f fc-log(l-l-z), k ^ 3.5, for z < 2. A low-luminosity 
cutoff was applied at A/|^™ = — 20 at z = for lo. 

^ In the general case [Sa^^f = 0.384^ [e^{Slof + [Sl^Y] with 
e from Eq.| and e = 1 - Ant/0.384. If e ~ 1 then Ant ~ 0. 
The same applies for aaox in El-M- 
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It should be noted that To and Tx are equivalent and 
interchangeable in this model. Eq. 6 can also be expressed 
in terms of Tx, in which case the X-ray luminosity function 
(XLF) has to be used instead of OLF. 

2.2. Monte Carlo analysis 

The Monte Carlo analysis was performed by generating a 
sample of quasars with z, Iq, and as follows: 
First, a redshift z was drawn from a given range (2:1,2:2) 
satisfying V{zi, z)/V{zi, Z2) = ri, where V{zi,z) and 
V{zi,Z2) are the volumes within {zi,z) and (2:1,22), re- 
spectively, and ri is a random number between and 1. 
The volume element in co-moving space is given by (for 

go = 0.5) 



dV{z) = 167r(c/i?o)^(l + z)"^-^(l + z- VlV^f dz. 



(8) 



Then, a second random number r2 was calculated and the 
intrinsic luminosity at the optical band Iq was determined 

such that /i'min 0o(L'o, 2) di'o / J-^^in 4>o{L' o, z) dL' o = 

r2, with the luminosity function (j)^{Lo,z) as specified 
above. After that, the expected X-ray luminosity Tx was 
obtained using Eq.|^ for a given Sox- The observed op- 
tical and X-ray luminosities, lo and 4c, are then drawn 
from the Gaussian distribution functions 5o('o — To) and 
ffx(^x — Tx), respectively. The corresponding fluxes were 
calculated from the luminosities assuming ao = —0.5 and 
ctx — —1.3 for the K-corrections, respectively. An object 
was accepted, if the fluxes were above the given flux lim- 
its. In general we simulated optically selected samples as- 
suming a threshold magnitude me^th = 20. If the objects 
turned out to have an X-ray flux below the assumed X-ray 
flux limit they were treated as 'non-detections'. 

The above procedure was repeated until a sample com- 
prising 300 objects was obtained. For each object, aox was 
calculated from the simulated lo and ic- We then investi- 
gated the ckox - relationship by means of a Spearman 
rank correlation test and least-square linear regression 
analysis. For a specific set of parameters of the quasar 
population, we carried out 10 independent trials and ac- 
quired 10 samples. All statistical quantities derived below 
are the mean of the individual values of the 10 samples, 
along with the statistical uncertainty for the mean. 

2.3. Results 

As the first and simplest case we assumed the intrinsic aox 
and (Tq;q^ in Eqs.||and|^ to be constant, i.e. independent 
of luminosity. In accordance with observational results we 
used clox = 1-40 and caox — 0.18. We considered 20 values 
of the i?cr-parameter from i?o-= 0.1 to 10, which are sam- 
pled evenly on a logarithmic scale. Eq.^ then determines 
the dispersions CTo and a^. We investigated optically se- 
lected samples with a size of iV = 300 and a limiting 
magnitude ms,th = 20. In a first attempt, we restrict our 



analysis to samples at a fixed redshift (z = 1), for which 
the results give an insight into the physical problem, and 
then we present cases for samples in larger redshift ranges. 

2.3.1. Sample at fixed redshift 

The Spearman correlation coefficients psp are plotted in 
Fig. 1 as open squares for 20 different Ra values, with 
the typical Icr error indicated for illustration. The prob- 
ability levels Pr — 0.01 and 0.001, at which the 'no cor- 
relation' hypothesis is ruled out, are shown at the cor- 
responding psp {N = 300) as dotted and dashed lines, 
respectively. Despite of the luminosity independence of 
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Fig. 1. Spearman rank correlation coefficients psp as a mea- 
sure of the Qox - lo correlation for simulated samples for various 
i?o-parameters; squares: sample at fixed redshift; filled circles: 
sample with redshift 0.2 < 2; < 3; open circles: sample with 
redshift 0.2 < z < 3 and containing only X-ray detections. 
The dotted and dashed lines indicate the pap for which the 'no 
correlation' hypothesis can be ruled out at the corresponding 
probability levels (one tail) of 0.01 and 0.001, respectively. 



the intrinsic mean aox the results show the emergence 
of a positive correlation with increasing T?^, which be- 
comes significant for iJ^- > 1. In Fig. 2 we show an exam- 
ple of an aox versus lo plot for a simulated sample with 

= 5. The dashed line represents a constant X-ray lumi- 
nosity (4c = 27.2 ergs~^ Hz~^). Obviously, the simulation 
indicates a much sharper cutoff in X-ray luminosity com- 
pared to optical luminosity, which has already been noted 
by Brinkmann et al. (1997) for observed data. 

Now we consider the effects of thresholds in the X-ray 
observations. At a redshift of z—1 a luminosity thresh- 
old of l^^ - 26.5 erg s^^Hz^i at 2keV indicated by the 
dotted line in Fig. 2 corresponds to a flux limit of ~ 
4 X lO^^'' ergcm^^ s^^ in the 0.5-2 keV band. Objects be- 
low this line have Ix > and can thus be 'detected' in 
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Fig. 2. ciox versus lo for a simulated quasar sample a.t z = 1.0, 
assuming a constant a^x and = 5. The lines indicate con- 
stant X-ray luminosities, 1^= 27.2 erg s~^ Hz~^ (dashed) and 
;x= 26.5 erg s^^Hz"^ (dotted). 

X-rays; whereas the rest with < l^^ are 'nopj-detections'. 
When only the 'detections' are consideredH, the X-ray 
thresholds can enhance the apparent a^^ - Iq correlation 
significantly, as can be inferred from the figure. 

The results are independent of redshift, which changes 
only the range of optical luminosities. 

2.3.2. Samples in redshift ranges 

We now consider samples with a redshift range of 0.2 < 
z < 3.0. The results of the correlation analysis are added 
to Fig. 1 for the case without (filled circles) and with an 
X-ray threshold (open circles). We find almost the same 
results as for the sample at fixed redshift. 

The slopes /3 obtained by fitting a linear aox f3 ■ lo 
relation to the simulated samples, are plotted in Fig. 3 for 
various values, for the case without (filled circles) and 
with (open circles) X-ray thresholds, respectively. The re- 
sults are fully consistent with those of the rank correlation 
test. Another effect to be noted is that the resulting av- 
erage Qfox of the obtained sample is consistent with the 
intrinsic aox {— 1-4) only for small (<SC 1), and it 
increases towards higher with a typical difference of 
- 0.1 for i?„> 1. 

The fitted slopes are found to be independent of the 
specified value of aox and insensitive to the X-ray observa- 
tional thresholds as well as to sample completeness. How- 
ever, they depend slightly on the dispersion craox (within 
plausible limits on cTao-^) in a manner that, at large Ra-, 
the slope increases with craox ^'^^ ^ given R„. We plot in 

^ Hereafter, we refer to such a case as 'with an X-ray thresh- 
old', otherwise we refer to data 'without an X-ray threshold' 



Fig. 3. Fitted slopes of Oox ~ /31ogLo for various Re values. 
Filled circles: no X-ray threshold considered and ctqox = 0.18; 
open circles: with X-ray threshold and ctqox = 0.18; triangles: 
with X-ray threshold and aaox = 0.15; squares: with X-ray 
threshold and craox = 0.20. The two shaded areas indicate the 
typical range of slopes (1 cr) given in two previous studies (see 
text). 

Fig. 3 the results for samples assuming Caox = 0-15 (trian- 
gles) and Ua^^ =0.2 (squares), respectively, ignoring the 
effect of X-ray thresholds. The shaded areas in the figure 
indicate the la confidence region of (3 obtained in previous 
studies {Pohs = 0.11 ±0.02, Wilkes et al. 1994, Yuan et al. 
1998; Pohs = 0.08 ±0.02, Green et al. 1995). It shows that 
the observed slopes can be reproduced within a relatively 
large region in the parameter space of Ra and cruox ■ fur- 
thermore, the slopes also depend only marginally on the 
low-luminosity cutoff Afg^'" chosen for To, which is some- 
what uncertain. 

In Fig. 4 we show the simulated aox and lo values for 

— 5. For 'non-detections' lower limits on aox are shown 
as crosses. This distribution, as well as that for samples 
at fixed redshift (Fig. 2), differs from what is expected for 
a linear aox ~ ^o relation, showing asymmetric structures 
with respect to a regression line. 

To test the correctness of the simulation procedure 
each sample was subjected to a self-consistency check us- 
ing the evolution weighted Schmidt's (1968) V/Vi-a, as in 
Chanan (1983). The simulated samples represented the 
envisaged luminosity function evolution as the resulting 
(V/Vni) were distributed randomly around 0.5 within the 
1 a uncertainties. 

We also performed similar analyses starting from Ix 
instead of To- After having determined z, Tx was drawn 
from the LF for Ix, 4'x(Jxj^)j which is approximated by 
the observed XLF. We used the broken power law LF as 
given in Boyle (1994), with a faint-end slope 71 = —1.6, a 
bright-end slope 72 — —3.3, and evolution rate k = 3.3 for 
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Fig. 4. ftox versus log Lo for a simulated quasar sample in the 
redshift range 0.2 < 2; < 3, assuming Ra = 5. Crosses indicate 
lower limits for objects which would not have been detected in 
the presence of an observational X-ray threshold. Thick bars in- 
dicate the best estimates of the mean ctox in lo bins which were 
obtained by taking into account non-detections (see Sect. 3.2). 

z < 2. Then, a corresponding To was obtained from Eq.[l] 
for the given clox- The observed luminosities were drawn 
in the same way as above. Similar results were found, i.e. 
an apparent aox - to correlation towards high Ra- . This is 
expected because of the striking similarity between the 
OLF and XLF (Boyle 1994). 

We conclude that an apparent oiox - correlation can 
emerge even for a population with an intrinsically constant 
cfox in the presence of large dispersions in the SED and 
for Ra> I- The degree of correlation depends on the i?^- 
parameter. Similar results hold for X-ray selected quasar 
samples, as well as for incomplete samples. 

3. Discussion 

3.1. Intuitive consideration 

The simulation results can be understood in terms of sim- 
ple intuitive arguments. We show schematic sketches for 
the Ix-lo relation in Fig. 5a and the corresponding aox - lo 
relation in Fig. 5b for a sample at fixed redshift. Following 
Sect. 2.1, the distribution of luminosities in the optical and 
X-ray regime is determined by the mean cfox- A constant 
clox = 1.4 is assumed (thick lines in Fig. 5a and b), as well 
as constant luminosity dispersions CTo and (indicated 
for an object at Q). The expected luminosities lo and Tx 
are distributed along the line of constant clox according to 
their distribution function (j)^(lo(Jyi), z), with a cutoff at 
the bright end of To and Ix for a sample with finite size. 

We consider the distribution of objects deviating from 
To and Tx due to luminosity scatter in the lx_-lo and aox - lo 
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Fig. 5. Schematic sketches for the l^-la and Qox-^o relation- 
ship. See text for a detailed explanation of the various symbols. 



planes. The 90% probability region for an object is con- 
fined within the two thin lines parallel to the mean aox- 
We consider two extreme cases: cto^CTx {Ra^ 1) and 
CTo^CcTx (-R(T <C 1). For CTo^CTx, the luminosity scatter is 
predominant in the optical luminosity, i.e. along the op- 
tical axis. The distribution of objects in the ^ - ^o plane 
forms a confined region below the horizontal dashed line 
AB, which corresponds to the line AB with slope 0.384 in 
Fig. 5b. An apparent correlation appears between aox and 
lo for such distributions, despite the intrinsically uncorre- 
lated relation. This effect explains the simulated aox ~ 
distribution in Fig. 2. On the other hand, in the case of 
CTo^cTx, objects are distributed in a region confined by 
line CD (dashed-dotted) and no aox - correlation is ex- 
pected from the high-Zo end. The actual degree of corre- 
lation increases with CTo relative to cTx, which explains the 
dependence aox-'o correlation on R„, as seen in Figs. 1 
and 3. Moreover, flux limits in X-rays (indicated by the 
line GH) , tend to enhance the apparent aox ~ correla- 
tion, as can be seen from Fig. 5b and Fig. 2. For an op- 
tically selected sample, the cutoff due to the optical flux 
limit (indicated by EF) systematically excludes objects 
with smaller aox values, which qualitatively explains the 
resulting increase of average aox for large R^. 
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We note that, in order to get a sharply confined aox - lo 
distribution at the high luminosity cutoffs, a steep lumi- 
nosity function for To (Jx) is required towards the high 
luminosity end. Our simulations show that the observed 
optical and X-ray quasar luminosity functions satisfy this 
condition. 

Using simulations basically similar to ours and a 
much steeper luminosity function (Gaussian distribution) , 
Chanan (1983) found the existence of a ofox ~ t correlation 
in his cases B and C, but not in case A, corresponding to 
1-4, Ra= 1, and i?^ ~ 0.8 in our analysis. These 
values are close to the critical point i?^ = 1, and the re- 
sults of Chanan are thus not a consequence of the reversal 
of dependent and independent variables in the regression 
analysis. Furthermore, the use of the i?o--parameter in our 
work, i.e. the relative strength of <7o and Ux, is physically 
more meaningful than the assumption that the luminos- 
ity of one wave band directly determines that in the other 
band (as in Chanan 1983). The case of Ra ^ 1 leads to a 
more pronounced high-luminosity cutoff for 4c compared 
to lo (see also Fig. 2). This effect is also seen in observa- 
tional data (Brinkmann et al. 1997). 

It can also be seen in Fig. 5a that, when performing a 
linear regression analysis for the Ix-lo relation, different 
distributions of the data near the highest luminosities can 
result in different slopes for different R^- Thus, a simple 
linear regression analysis method is not always adequate 
to quantify a luminosity correlation, especially for data 
with large inherent scatter. 

3.2. Dependence of on — intrinsic or apparent? 

Our results show that an apparent dependence of ctox on 
lo can emerge from data for a quasar sample with no in- 
trinsic dependence. Thus, the explanation of the observed 
CKox lo correlation as a physical relation in quasars, as 
taken for granted in previous work, must be questioned. 
Although the current study does not allow to unambigu- 
ously distinguish between an intrinsic or an apparent de- 
pendence of Qfox on lo, we may reach some conclusions by 
simple considerations. 

If our model is a good description of the quasar lumi- 
nosity correlation and dispersion, we expect the existence 
of an apparent correlation to some extent, unless in cases 
of rather small R^ 0.3), i.e. when the dispersion in lo is 
at least 3 times less than that in 4c, which seems unlikely 
considering the diversity of the big-blue-bump component 
in quasar energy spectra. In fact, the aforementioned char- 
acteristics of a much tighter high-luminosity cutoff for 4c 
than for lo seems to exist in most quasar samples of con- 
siderable size. Thus, an intrinsic aox ~ lo dependence, if it 
does exist, must be weaker than it appears from the data. 

To quantify these effects, we carried out similar simula- 
tions incorporating an intrinsic luminosity dependent Sox 
as in Eq.^, /Sjnt ^ 0, and compare the obtained slopes /3 of 
the aox - lo relation with the observed values. For each grid 



point in the Pint— Ra parameter space, we repeated the 
simulations for 200 trials and counted the number m of tri- 
als for which (3 fell into a region around the observed value 
/3obs ± ujcr, at a confidence level ^. The overall confidence 
level of reproducing the observed slope is thus ^ • m/200, 
if systematic uncertainties, such as sample incompleteness 
etc., are ignored. For Pohs — 0.11 ± 0.02 (Wilkes et al. 
1994, Yuan et al. 1998), Fig. 6 shows the contours of the 
68% (thick lines) and 90% (thin lines) confidence regions 
of the joint Ant — Ra distribution for Uaox =0.15 (dotted), 
0.18 (solid), and 0.20 (dashed), respectively. The results 
show that only a weak dependence {Pint ^ 0.05) is needed 
for <; 1, when the dispersion in lo is comparable to 
or larger than that in 4c; and almost no intrinsic depen- 
dence is needed for larger i?^ 3. It also shows that the 
most probable (68%) region appears at high R„ (^ 1), low 
/3i„t(^ 0.07) values. 
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Fig. 6. Contours of the 68% (thick lines) and 90% (thin 
lines) confidence region in Re and /3int, the slope of the intrin- 
sic Qfox-^o dependence, for auox ~0-15 (dotted), (JQq^=0.18 
(solid) and craox =0.20 (dashed). 



More qualitative evidence comes from some features 
of the observed aox - relation. Recent studies of large 
quasar samples (Avni et al. 1995, Brinkmann et al. 1997, 
Yuan et al. 1998) have shown that the observed aox-lo 
correlation shows a more complex behavior, with only 
weak or even vanishing correlation at low optical lumi- 
nosities. We find that such features are expected, as the 
apparent correlation for R„ > 1 arises mainly from objects 
with high optical luminosities. This property has been ver- 
ified by estimating the average aox in different lo bins for 
the simulated data. This can be seen in Fig. 4, where we 
plot the average aox (thick bars) for five lo bins, which 
were obtained by incorporating upper limits for the 'non 
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detections' using the maximum likelihood method devel- 
oped by Avni et al. (1980). 

A comparison of the observed OLF and XLF and their 
evolution may also give, in principle, constraints on the 
intrinsic cfox-^o, or Ix-Io relation (Eqs.^ and |^). For a 
power law luminosity function for Lq with index 70, and 
luminosity evolution ^ (l-|-z)'^°, it can be shown (see also 
Kriss & Canizares 1985) that the corresponding LF for Lx 
is also a power law (index 7x) with evolution (1 + z)''^ and 

7x - 1 = (70 - l)/e, /cx = e/co. (9) 

For a Gaussian approximation of the luminosity disper- 
sion, the luminosity functions 4>a(lo, z) , (/'x(^x, z) have the 
same forms and evolution as the observed luminosity func- 
tions (j)o{lo, z), 4>x{lx, z). Thus, e and /3int can be estimated 
from the observed ko, fcx, 7o and 7x. 

The optical evolution index is found to be fco — 3.45 ± 
0.1 (Boyle 1994). However, the X-ray evolution index is 
somewhat uncertain — ranging from fcx — 2.56±0.17 (Delia 
Ceca et al. 1992, qq = 0, Einstein EMSS data only), fcx = 
S.Otoa (Jones et al. 1997) and fcx = 3.34 ± 0.1 (Boyle et 
al. 1994) to fcx = 3.2 - 3.5 (Franceschini et al. 1994). The 
latter three values were obtained by incorporating ROSAT 
data, and they are consistent with e = 1 and /3int — 
within their 2 a errors, i.e. no intrinsic dependence of aox 
on lo is needed. Even if we take the smallest value of fcx = 
3.0 of Jones et al. (1997), we have e = 0.87 and /3int = 0.05, 
implying only a marginal dependence, weaker than the 
previous claims of (3 ~ 0.1. The bright-end slopes of the 
XLF (72 = 3.30l^;i^g) and of the OLF (72 = 3.9 ± 0.1, 
Boyle 1994) are consistent with each other within their 
2(7 errors, and thus also consistent with e = 1 and /Sjnt = 
0. We conclude that a constant cfox independent of lo is 
not ruled out by a comparison of the observed luminosity 
functions in the optical and the X-ray band. 

4. Conclusions 

We have performed Monte Carlo simulations to study 
the luminosity correlation between the optical and X-ray 
bands for quasars. We have used a generalized model in 
which the luminosities in the two wave bands are rep- 
resented in terms of the respective expected luminosities 
with dispersions {lo = lo + Slo, l^ — Ix + Six)- We have 
shown that the increase of aox with Lo (equivalent to 
Lx oc Lq with e < 1), as found in observational data, 
can emerge in a sample with an intrinsic luminosity inde- 
pendent Qfox (or Lx oc Lo), provided that the dispersion of 
the optical luminosities deviating from the average SED 
are similar to or larger than that of the X-ray luminosi- 
ties. Our simulations verified the results of Chanan (1983), 
which were achieved for special assumption about the lu- 
minosity functions. We suggest that the observed Uox-Lo 
correlation is, at least to a large extent, apparent and not 
necessarily an intrinsic property of the quasar population. 



Our model is more general than previous considera- 
tions. For ^ 1 (implying lo ~ lo and (Jaox ^ 0.384crx) 
the model reduces to the commonly used description, in 
which ^o is the primary luminosity and the dispersion in 
the SED is attributed to that in the X-ray luminosity. The 
same holds for ^ 1, but with interchanged roles of lo 
and Ix- We argue that the effect of the relative strength 
of the individual luminosity dispersions in the two bands 
should be taken into account in analyses of quasar lumi- 
nosity correlations. Since the arguments are valid for any 
other two wave bands, we expect this effect to play a role 
in luminosity correlations between radio, infrared, opti- 
cal, and X-ray wave bands as well. The determination of 
the _Ro.-parameter, and thus of the luminosity scatter in 
the individual wave bands, is important to understand the 
broad band emission of quasars. 

We finally note that a similar effect as presented for the 
Q!ox — correlation might also apply for the well-known 
Baldwin effect, i.e. the inverse correlation of optical emis- 
sion line equivalent width with optical luminosity (Bald- 
win 1977). Since the equivalent width basically can be re- 
garded as the ratio of two luminosities (emission line and 
underlying continuum), the structure of the problem is 
similar to the one presented in this paper. This is particu- 
larly interesting, because there still is no accepted physical 
explanation for the Baldwin effect. 
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